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ABSTRACT 

1 2 

Cholshevnikov 5 has published the results of analytic studies demonstra- 
ting that bounds on the zonal geopotential coefficients decrease at least as rapidly 
as l/n and bounds on the tesseral coefficients at least as rapidly as l/mn . 

We show that these bounds are almost certainly conservative, and then proceed 
to develop a number of possible modifications designed to incorporate extensive 
data on the Earth's density distribution into the analysis. None of the modifica- 
tions are implemented; the most promising is, however, discussed at some 
length. 
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I. 


INTRODUCTION 


The normalized geopotential coefficients can be given in terms of the density 
of the Earth, p( r , z ,<p), and the normalized associated Legendre functions, P (z), 
used by geodesists, as 




nm 


, „1 2tr A(z,cp)_ n+2 

P nm <2) {r^}] d f ,ds dr 
-10 0 A ,1 

e (1.1) 


where 


A 

e 

z 


radius of the smallest sphere centered at the origin 
and circumscribing the Earth 

cos 0 ; 0 = polar angle of spherical coordinates 


<p = longitude 

A(z ,cp) = distance of that point on the surface of the Earth, 
defined by z and <p , from the origin 


Knowledge of the behavior of these coefficients as functions of n and m would be 
useful in many ways, particularly in providing a rational basis for truncation of the 
geopotential in applications to orbital mechanics, and in the formulation of proce- 
dures for geopotential modeling. 


Since the density, p, is an empirical quantity, an analytic evaluation of 
the coefficients from Eq. (1.1) is not possible. With sufficiently detailed density 
data, numerical integration of the right hand side of Eq. (1. 1) could be used to con- 
struct a table of values for the coefficients. Although sufficient data for such a 
calculation are not presently available, a considerable amount of data does exist, 
on the basis of which various density models for the Earth have been constructed. 
The purpose of this report is to discuss how the available data and models might 
be used to investigate the behavior, not of the coefficients themselves, but of bounds 
on the coefficients. It has been generally supposed for a long time that the geo- 
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potential coefficients decrease as n increases, a conjecture supported by em- 
pirical estimates based on examination of tables of empirically determined coef- 
ficients. The most sophisticated analytic estimates, to my knowledge, are those 
12 

of Cholshevnikov, ’ who gives four mathematically rigorous derivations starting 
from Eq. (1. 1). His results are on the rate of decrease of upper bounds on the 
magnitudes of the coefficients, and we shall describe them in some detail in the 
next section. Here we wish only to comment briefly on the approach, and intro- 
duce the mathematical background necessary for his development. 

2 

A bound of, say, A/n for I C I is not the same as a decrease of I C I 

2 nm 1 1 nm 1 

as A/n ; |c I might decrease more rapidly; or, perhaps, starting from a 

value much lower than the bound, it might increase for awhile; or it could even 

have a damped oscillatory character. However, whatever the variation of I C I 

2 nm 

with n , I C 1 would be constrained to lie below A/n , and the bound would 
nm ' 

imply lim C = 0 . 
n_z nm 

Cholshevnikov' s bounds are obtained by using standard inequality relation- 
ships on the absolute value of an integral, bounds on P (z) and/or its integral 

nm 

from -1 to +1 , and, for the sharper results, a mean value theorem for integrals 
based on the concept of a monotone function. A monotonically increasing function 
of x satisfies the inequality 

f(Xj)*f(x 2 ) for x 1 >x 2 (1.2) 

It differs from a strictly increasing function since there may exist one or more 
intervals over which it remains constant. A monotonically decreasing function 
is defined by reversing the second inequality. We now state the mean value theo-- 
rem; a proof will be found in Reference 3. 
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Theorem: 


If F(x) and G(x) are continuous in the interval 
(a,b), and also if F(x) is monotone (either in- 
creasing or decreasing) in this interval, then 
there exists an x^ in the interval such that 


*1 


J F(x) G(x) dx = F(a) J G(x) dx + F(b) J G(x) dx 


(1.3) 


Clearly, the bounds on C and S must involve the density function 

nm nm 

and the shape of the Earth defined by A(z,c p), another empirical quantity. Chol- 
shevnikov eliminates A(z,<p) by using A g as the upper bound for the integration 
with respect to r in Eq. (1.1). This implies that the density function must be 
set to zero in the gap between the surface of the Earth and its circumscribing 
sphere. Cholshevnikov incorporates the density function in his analysis through 
the "global' r functions of density 


p , a global constant (p . = 0) 
Tnax "min ’ 




r . 5p . * 

v^(r, 0, <p) = J | — | dcp = variation of p with respect to (p (1.4) 


0 , 
.z 


v_(r, z. <p) = J |-^|dz = variation of p with respect to z 


Clearly v^ and y^ are monotonically increasing functions of <p and z , respec- 
tively, and a little thought will show that (v - p) and (v - p) are also monotonically 

(p Z 

increasing with respect to <p and z , respectively. These comments perhaps al- 
ready suggest how Cholshevnikov derives his bounds. Since several modifications of 
his analysis are almost immediately apparent, we shall devote Section n to an ab- 
breviated derivation of his results which, while rigorous, have two defects: 


* For a more rigorous definition of "variation," see Reference 4. 
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1. One would conjecture that they must be fairly conservative 
since the fact that the integrand in Eq. (1. 1) has one sign 
over approximately half of the region of integration, and 
the opposite sign over the rest, appears to have been fully 
exploited only in his Theorem II. 

2, The parameter, n , can be reasonably well estimated. 

'max 

Adequate data for the estimation of bounds on v^(r, 0, 217) 

and v (r, 1, <p), used in the sharper results, are probably 
z 

not available. 

In Section III we discuss some possible modifications of Cholshevnikov’s 
methods, aimed at eliminating the defects mentioned above. The primary change 
is to introduce- local formulations rather than the global formulations implemented 
by Cholshevnikov. In Section IV, we select what appears to be the most promising 
of the local formulations for a more detailed analysis. It is beyond the scope of 
this report to implement any of the suggested formulations. We do indicate, how- 
ever, how data on the density function and the shape of the Earth might be used in 
an implementation. In the last section, we make a few comments on "worst case" 
constructions and their role in deriving bounds on the geopotential coefficients. 
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II. 


CHOLSHEVNIKOV’S BOUNDS FOR GEOPOTENTIAL COEFFICIENTS 


The first of Cholshevnikov's papers gives two bounds for the coefficients 
of the zonal harmonics, one in terms of P max and the other (which is sharper) 
in terms of the variation v . The second paper carries through the corresponding 

Z 

analysis for the coefficients of the tesseral harmonics, the major difference in the 
two investigations being that the variation v^ is used in the second. We shall state 
these four theorems and sketch the proofs, modified to apply to normalized coef- 
ficients (Cholshevnikov used conventional associated Legendre functions and un- 
normalized coefficients). The analysis applies to the defining equation (1. 1), for 
the geopotential coefficients, written in the form 


e n+2 
r 


.1 .2 it 


vu.i 


nm 


-1 0 


This form differs from (1. 1) in the upper limit of the r-integral. There is no loss 
of generality; the replacement of A(z, <p) by A can be compensated by setting 
the density p to zero in the gap between the actual surface of the Earth, defined by 

r = A(z,<p) (2. 2) 


and the circumscribing sphere of radius A , since from the definition of A 

6 Q 


A( z , <p) ^ A 


(2.3) 


Theorem I: 


4tt p A 
I ^ 'max e 

nO 1 (2n+l)(n+3) 


(2.4) 


i 
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Proof: 


From Eq. (2. 1), S n ^= 0 and 


e n+2 , 1 2 n 


C £ 
nO 2n+l 


SlJ J lp( r > Z, <« H P nln <z>l d 0 <fc} <ir 


-1 0 


2 it p _ • n+3 As „1 

Jl^l- 

(n+3)A g 0 -1 


from which the result follows immediately on use of the inequality 

.1 

l a0' 


I ' P i0< z > 


dz ^ 2 


which Cholshevnikov credits to Hobson. Note that the second inequality 
in (2. 5) is strictly "less than" because if p = P max for all r , z , <p, 
Eq. (2. 1) implies C = 0 for all n except for n = 0 . 


4 

4 7T A 3 

Theorem II: I C I =£ 

1 nil * 


2 v + p 
z max 

max 


“0 < 2n+1 > < n+3 > (2n+l) 


where 


v = LUB(v (r, 1, c p) for all r, <p) 
z z 

max 

3 

Note: This theorem implies that LUB|C J decreases as const/n 

for large n. Cholshevnikov shows, in his first paper, that this 
theorem is a "best" result in the sense that a density distribution can 
be constructed for which |c itself decreases as l/n 3 . We com- 
ment-inmore detail onjthis in:Sectiqn. V. 


Proof: 


Rewrite Eq. (2. 1) in the form 


(2.5) 


(2.6) 


(2.7) 


( 2 . 8 ) 
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1 r e r n+2 r 2 ,I rr 1 

C nO = iST J ~iT J {J , v z< r ’ z ' P nO (z > dz 

0 A 0 -1 

e 


- J [v_(r, z, <p) - p(r, z,<p)]P n0 (z) dz}d<p dr 


Since both v and v - p are monotonically increasing functions of 
z z 

z , the expression in "curly" brackets may be replaced by 

r 1 r 2 

v z ( r> 1,<p )J P n0 (z) dz " ^-P( r >" 1 » p n0 ( z ) dz 

z. -1 


- [v z (r, l,<p)-p(r, l,<p)] J P n (z) dz 


(2.9) 


( 2 . 10 ) 


where use has been made of the mean value theorem, stated in the 
Introduction, and of the fact that 


v,(r, -1, <p) = 0 


( 2 . 11 ) 


Further, since 
.1 


J P n (z) dz = 0 = J P <z) dz + J P n (z) 

-1 -1 z„ 


dz 


( 2 . 12 ) 


we can combine the last two integrals of Eq. (2. 10) to obtain 

v r -- 1 -*»[J l P n (z) dz " J p n ( z ) dz ] 

Z 1 Z 2 

+ [p (r, 1, <p) - p (r, -1, <p)] J -P-(z) dz 


(2. 13) 


for the expression in curly brackets. Cholshevnikov now makes use 
of another inequality 
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I J P n (z) dz j £ 


4 

2/22 


V 7T (n-1) (2n+l) 


(2.14) 


to obtain 


2 TT 


_ i^-J- r i 

nO 1 2n+l 


2^22 


2 v + p ] - 

z max / 

max V tt (n-1) (2n+l) 


d <p dr 


(2.15) 


from which the desired result follows immediately. The reader is 
referred to Cholshevnikov's paper for the derivation of the inequality 
(2.14), Note that we have dropped the subscript 0 on P , since 
for m = 0 , the associated functions become Legendre polynomials. 


Theorem III: 


f run' *1 

Us | J 

nm 


An J 2 p A 3 
max e 

(2n+l) (m+3) 


for m>0 


(2. 16) 


Proof: 


From Eq. (2. 1) 


e n+2 1 2 tt . 

C + i S = —77 P (z) e im< ^p(r, z, <p)d<p dz dr 

nm nm 2n+l „n J, nm' ' LJ„ r/ r J 

0 A e - 1 0 (2.17) 


and hence 

lc 


,\C L • , „e n+2 „1 _ „2 tt 

{ |s nm |} S |c nm + 1 S nm> S i^T I V I, P n m< z) [I Z ’ <P) M dz dr 

1 nm 1 0 A -1 0 


2 tt p _ n+3 _ e 

'max f r 


(2.18) 


2n+l 


[_J_] 

1^ -J A 


(n+3) A^ 0 


where use has been made of another inequality credited to Hobson: 

.1 

|] 

nm' 


p-L 

J l P mn (Z >l dz S 


(2. 19) 
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Theorem IV: 


r 'Wi 

1 Is li 


4*/2 v A 
0 e 
max 


nm 


m(n+3) (2n+l) 


( 2 . 20 ) 


Proof: 

lc 


e n+2 .1 


2 77 


(|s nm |} ^I 0 Vj l^lf P<r -^ , {“n“^ d " |dZdr 

1 nm 1 0 A -1 0 


( 2 . 21 ) 


This time, using the variation with respect to cp, we obtain a bound 
on the (p integral: 


| f 27T p{ C ° S m<P }dcp | = |v (r, z, 277) n C ° S m ^)dcp 
^lsinm<pJ 1 <p' tsinmcpJ ^ 


,2 7T 


‘Pi 


(2. 22) 


+ P(r. 2 . 0 )J '{“ S ”^}d?>- [v(r,z > 2 ff )-p(r,z,2 1 r)]J {“ S ”,£}<M 
0 *2 


where we have used v^r, z, 0) = 0 . From continuity, we have 


p(r, z, 0) = p(r, z, 2 77) 


(2. 23) 


and hence 
,27 7 
0 


<Po 


2 77 , 


■I p {stam£} d '* ,| = >V r ’ Z ’ 2,,> I {stom^} d?,+ P(r ' z ' 0) I 


<p. 


2 v 


(2. 24) 


v (r, z, 2 77 ) , sin mcp., - sin m<p, ^ 

0 K ’ f ‘2 1 1 j ^ max 


<P„ 


f 2 ^11 . 

m \cos mcp - cos mcp J m 

1 Lt 


Proceeding now as in Theorem El, the desired result is easily obtained. 


We now combine the methods used by Cholshevnikov for Theorems n and m 
to obtain a fifth theorem: 
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Theorem V: 


J nm \ 

l ie I J 


nm 


2 ir A 

e — [ 2 v 

(2n+l) (n+3) z 


+ A 


max 


max 


Hi 

-1 


P (z) dz 
nm' 7 


(2. 25) 


Proof: 


e n+2 „2 tt 


l°«n + * S nml S 2^1 I V I (e 1 J p(r. z, <p) P nm (z) dz| d<p dr (2. 26) 

0 A 0 -1 


or 


e n+2 „2ff 


{ le |} " n~7T [ [ £2 v + p ] • P (z) dz dep dr 

t. S J 2n+l J. A n J. z 'max 'J, nm' ' 1 ^ 

1 nm 1 0 A 0 max -1 


where the z-integral has been treated as in Theorem n. Note that use 
of inequality (2. 20) would yield a result inferior to Theorem IV. Pre- 
sumably, Cholshevnikov would have included this theorem also had an 
estimate of ilj p nm ( z ) dz | decreasing faster than l/m been available 

to him. Analytic bounds on P (z) dz are not known to me. 

il nm' ’ 


Of his various theorems, Cholshevnikov proved only that Theorem n is a 

best result, i. e. , that density functions (however "unrealistic") exist for which the 

equality, rather than the inequality, holds. If Theorem n is "best," Theorem I 

cannot be "best." In fact, since both Theorems I and in involve only p and 

'max 

bounds on J l p nm ( z ) I dz (which do not exploit the oscillatory character of p nm )> 
both can certainly be sharpened. Theorems IV and V make use of a variation of 
p and include the oscillatory character of both P^ (z) and cos (mep) or sin (m<p); 
it is possible that at least one of them is a "best" result. The variations v Zmax 
and must, however, be enormous numbers, and even rough estimates of 

their values are probably impossible to make from available density data. In 
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spite of this, Theorem II, involving v„ is a "best" result, at least for the 

z max 

order of magnitude of the rate of decrease. We shall return to this point in 
Section V. 

The theorems discussed in this section all involve constant factors de- 
pending upon "global" functions of the density: p , v„ , and v ^ . In 

the following section, we propose a method for sharpening at least the constant 
factors by constructing a local rather than a global formulation, with the possi- 
bility of using local rather than global bounds on the density and its variations. 
Such a formulation should enable us to make more effective use of the available 
density data. 
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m. 


LOCAL PROCEDURES FOR BOUNDING THE GEOPOTENTIAL COEFFICIENTS 


Our local formulations will be based on rewriting the defining equation (1. 1) 
for the geopotential coefficients in the equivalent form, for m> 0 : 

C n-m+1 2m-l 

{ s nm hiM .1 <*•« 


nm 


■ 1=1 3=1 


with 


where 


Z i+1 ^j+1 


I.. 

U 


-/ J 


_ , , f cos mtp'i „ , 

P (z)j . ■ y R (z , <P) dip dz 

' ' l sin m<p J n ■ ^ 


*i 


nm 


(3.2) 


„A(z , <p) 


n+2 


» r / -£ — 

R(z,<p) = p(r , z ,<p) — — dr , 

n J 0 A n 

e 

the z. are the n-m+2 zeroes of P (z) with z , = -1, 
i „ nm' 1 

Z n-m+2 + ’ 

, _ r cos m cp 1 , , 

the <p. are the 2m zeroes of i . 1 in the interval 

j n ^ t sin mto J 

J 0 ^ <p< 2n . ^ 


(3.3) 


The case m = 0 requires special treatment and, since it may not be amenable 
to a local formulation, it will not be discussed in this report. 


The z. and <p^ define a grid of spherical quadrilaterals (triangles at the 

poles) bounded by the latitude and longitude lines corresponding to the zeroes of 

that spherical harmonic which defines C (or S ). Note that this representa- 

nm nm 

tion is unique for each of the coefficients C , S as n and m run from 
nm nm 

zero to whatever upper limit is imposed. 
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The particular representation of Eq. (1. 1) as a double sum has the property 


that for 


z. ^ z ^ z. , P does not change sign 
l i+l nm 

(p. ^ to ^ (D. , j °? S m< ^ f does not change sign 
J j+1 fsinmc pi 


(3.4) 


which implies that the inequality usually present in relating the absolute value of 
an integral to the integral of its absolute value degenerates for I.. to 


I.. = ± 
ij 1 


J I V I { j sin mp j} V z • ^ dz 

1. <p. I V- I 

l J 


(3.5) 


Further, the fact that the z. and </j. are zeroes of P (z) and { C ? S m< ^T . 

l j nm ' l sm m<p J 

respectively, implies that 


I.,, . and I. have signs opposite to that of I.. 
i±l,] i,j±l & w ij. 


(3.6) 


since one factor of the integrand changes sign if either i or j (but not both) 
change by 1 , while the other factor does not change sign. Similarly, 


I. , , . , , has the same sign as I.. 


(3.7) 


for all four combinations of the +'s and -'s. 


Now let us examine the contribution to the sum over i and j in Eq. (3, 1) 

of two neighboring sectors with common boundary . Referring to the sketch, 

we see that the positive contributions come from one zigzag path, and the negative 

contributions from another. Note that if R were constant, the total contribution 

n 

of this pair of sectors would vanish and, hence, C and S would also vanish 
* . nm nm 
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since there are 2m sectors and, therefore, an 
integral number m of pairs. Next let us set 


R.. = LUBR(z,ifl) 

i], max n ^ 


(3.8) 


R. . . = GLB R (z ,<p) 

ij,mm n 


3+1 


for "block i , j", defined by the inequalities of 
(3.4). : 
sketch. 


(3. 4). If we take I to be negative, as in the 


I. . .. + I. . < (R. . -R.. .)|"f P (z) dz][ f 3 { C0S m<p }d(pl 

i,j-l i,j i],max ij,min LJ^ nm JLJ_ Lsinmip J 

(3.9) 


<P* 

f 

<P 


3-1 


2[R.. -R.. .] 

VI , max ij , mm 


m 


P (z) dz 
J nm' ' 


z. 

l 


and the more nearly equal are R. . and R. . . , the more fully can the 

lj, max 13, mm 

cancellation effect from integrands of opposite signs come into play. Also, the 

more accurately the interval integrals between successive zeroes of can be 

estimated, the more accurately can the n , m dependence of the bounds of the 

geopotential coefficients be determined. References 5 and 6 discuss various 

properties of P nTr| ( z ) pertinent to this latter problem. Our primary concern 

here is with the incorporation of data into the analysis — that is, with possible 

procedures for the estimation of the bounds R.. and R.. . of R (z , <p) 

13, max 13, mm n 

for the i , j block. 


We first recall that R depends upon the two empirical quantities; A(z ,<p), 

n 

which defines the shape of the Earth, and the density function p(r,z ,<p ); and 


14 



that, in fact, these empirical quantities appear explicitly only in in this for- 
mulation. Basic inequalities involving A(z ,<p) which might be useful are 


A . < A(z ,<p) ^ B(z) ^ A ' (3.10) 

where B(z) is the distance from the origin to a point with polar angle arccos z 

on the smallest oblate ellipsoid with semi-major axis A g , centered at the origin, 

and circumscribing the Earth. A . is the radius of the largest sphere, centered 

mm 

at the origin, and contained by the surface of the Earth. The inequalities (3. 10) 
are "global." Local LUB's and GLB's of A(z,< p) would clearly be useful, 
and can probably be fairly well sized over most of the globe. Probably the real 
key to the estimation of realistic bounds for R^ lies in finding an effective way 
to handle p. A discussion of this problem is given in the next section for the 
"radial" formulation developed above. 


Before discussing the radial formulation in depth, however, we wish to 
point out that the local approach, including the symmetry properties of the "block" 
integrals incorporated in Eqs. (3. 5), (3. 6), and (3. 7), can also be given the fol- 
lowing formulations: 


with 


z. 


n+2 


i+1 B(z) 

I.. = | | P (z) — $.(r,z) dr dz 

ij J J n nm v A n 


A 

-e 


\ n+2 +1 


(3. 11) 


r r__ p {cosm^j 
. n J IsinmtpJ i 


0 A <0. 

e 


<P 


j+1 


Ojlr.z) = J P(r^.«’){3° n S “^}d'P 


(3. 12a) 


p i 
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(3. 12b) 


The inequalities (3. 10) have been used in writing Eq. (3. 11) and, in order to main- 
tain Eq. (1.1), p(r , z , <p) must be set equal to zero in the gap between the surface 
of the Earth defined by A(z , <p) and the surface defined by B(z) in the first of Eqs. 

(3. 11), and the circumscribing sphere of radius A in the second. These formu- 

© 

lations lend themselves to the application of the methods used by Cholshevnikov in 

a local rather than a global manner; that is, p and the variations v and v 

max <p z 

can be defined for the individual blocks to obtain bounds on $ and Z. . 

J i 


■; 


i+i 


z. 

1 


p(r,z,<p)P nm (z) 


dz 


Proceeding just as in the analysis leading to inequality (3. 9) and using the 
symmetry properties of the blocks, we find the corresponding inequalities 

r n r B(z) n+2 Z i+1 

I. . 1 +. I.. < - $.. . dr • P (z) dz 

i,j-l ij L ij,max vj,minJJ n J nm' 


z. 

i 


A 3 Z -,1 

A _ i+l 
e 


$. . - $. . .1 f I P (z) I dz 

L ij,max ij,mmJ n+3 J 1 nm 1 


(3. 13) 


z. 

i 


and 


e n+2 


I. . , + I..<fz.. - Z.. .If dr • 

i,]-l 1 3 L ij, max ij.mmJ A n 




j 

J | cos m<p 




sm m cp 


}d<p 


= [z.. - Z.. . 1 

L ij , max ij,minJ 


2A 

e 

m(n+3) 


(3. 14) 


The task of estimating bounds for Z. and <5^ appears to be somewhat less 
tractable than that of estimating bounds for . Further, in the first of these for- 
mulations, (3.13), the empirical function B(z) is present in the first of the inequalities 
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although, as indicated in the second inequality, it can be eliminated at the expense 
of weakening the bound. 

No effort has been made to implement these two formulations; they are 

suggested because the utilization of the variations and by Cholshevnikov 

appeared rather effective, although the difficulty of estimation of global bounds 

for v and v renders his results of theoretical rather than practical interest. 

<p z 

The estimation of local bounds, particularly for n»m»l (so that the blocks 
are all small), may be feasible. 



IV. 


INCORPORATION OF DENSITY DATA INTO A LOCAL RADIAL 


FORMULATION 


We designate our principal local formulation, based on R^ , as a local radial 

formulation for rather obvious reasons: it is local and R is the radial part of the 

n 

defining integral (1.1) for the geopotential coefficients. The function R^ contains 
all the empirical quantities and, once we can establish bounds on R^ for each block 
i , j , the bounds on the geopotential coefficients are constructed from a weighted 
sum of the interval integrals 



z. 

i 


P 

nm 


(z) dz 


(4.1) 


of P nrn ( z ) between successive pairs of its zeroes, as outlined in the previous 
section. A computer program 6 has been written to calculate the interval integrals 
and can operate with accuracy adequate for this purpose up to degree and order 
40 , 40 . Analytic estimates of the interval integrals as functions of n and m 
would be highly desirable. No such estimates are known to me, and their deriva- 
tion appears to be a problem of considerable difficulty. 


The local formulation (3. 1) of the defining equations for the geopotential 
coefficients has somewhat the character of a numerical integration. The original 
triple integral of Eq. (1. 1) has been replaced by a double sum of block integrals: 




nm 


n-m+1 2m-l Z {+\ ^i+l 

t IJ/ P nm (z) {s°nm^} R n (z ' <p>d,i,d!; < 4 ’ 2 > 


2n+l 


i=l j=l 


with 


A(z ,<p) 


n 7 r I 

R (2><P) = p(r , z ,<p) 

n 


n+2 


n 


dr 


(4.3) 


18 


No approximations have been introduced up to this point. Notice that the sizes of 

the blocks are determined by the values assigned to n and m , and may be too 

large for a conventional numerical integration. This could, of course, be remedied 

by further subdividing each block. The available density data and data on the shape 

of the Earth's surface defined by A(z , cp) are, however, not adequate to carry out 

a direct numerical integration of Eqs. (4. 2) and (4. 3) with sufficient accuracy to 

give realistic estimates of the values of C and S . We seek, instead, bounds 

nm nm 

on R to be used as input for calculation of bounds on the coefficients, 
n 


One way to attack this problem is to note that r is a monotonically in- 
sreasing function of r and hence, using the mean value theorem in the Introduction, 

n+2 „ r n^’^ n+2 A(z,<p) 

R n ( z ><P) = [“ — | p(r , z ,<p) dr + £ f 

A 0 


A 0 
e 


£ \ > 'T / — 

p(r,z,(p) dr 

/ rr tr\\ ~ 


y z »<p) 


.n+2 4 A(z,<f>) 

= a (z ,(p) r 

A n r (z , <p) 
e n ^ 


p(r,z,<p) dr 


(4.4) 


bis approach has the consequence of eliminating the factor — r which we are now 

n+o 

jcustomed to associate with the radial integral. This loss may, however, be eom- 
msated by the lower limit "r^z , <p ) on the last integral. If we transform the vari- 
>le of integration by 


e 


r (z , <p) 

< 4 - 5 > 

P A 


n R (z ,<p) becomes [Eqs. (4.3) and (4.4)] 

£„„Jz ,<p) 


>ax'"’ r ' ’max (z,<p) 

R(z,«?)= A 0 J p(4,-z,<p)4 n+ -d4-= A- p(r,z, <p) d<p 

0 £ (z.v) 


C n (z,<p) 


(4.6) 
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We note that one factor, at least, of the integrand of the first integral decreases 
rapidly with increasing n over most of its interval of integration. On the other 
hand, £ max » the corresponding factor for the second integral is very close to 
unity and will decrease comparatively slowly for values of n of any immediate 
concern: for example, k must exceed 128 before . 9966^ drops below \ — 

. 9966 is approximately the ratio of the polar to the equatorial radius of the Earth. 
One would thus expect 4 n ( z ><P) to approach £ max ( z »<°) fairly rapidly as n in- 
creases. To translate this expectation into a quantitative estimate of 4 n ( z ,<p) is 
the difficult step in calculating bounds on R^z.o). In principle, £^(z ,<p) could 
be evaluated by equating the two integrals in Eq. (4. 6) and solving for f^(z , ip) , 
but this is not feasible in practice. 


We postpone, briefly, the discussion of this problem in order to investigate 
just what information on £^(z , <p) might be required in order to place bounds on 
R q (z ,<p). Using the second form of R^(z , <p ) given in Eq. (4. 6), we see that if 


A (z , <p) = £ (z , <p) - 4 (z , <p) 
n max n 


(4. 7) 


is "small, " fairly tight bounds can be placed on R^(z ,<p) as follows: 


A e Cax^ Z ’ •P) P mia (z.V)<\(z.‘P)<\ 

Cax^ Z ’ ^ ^n^ Z ’ ^ ^rnax^ Z ’ ^ 


(4. 8) 

where p (z ,<p) and p . (z , <o) are maximum and minimum values of the density 
'max' 'min' _ 

along the " ray" (z ,<p) in the shell bounded by 4 (z , <p) and 4 (z , <p). Now re- 

ULRX U 

call that what we need for Eq. (3. 9) are the bounds R.. and R.. . on 

' ' ij , max 13, mm 

R n (z,<p) for each pair of blocks i , ( j— 1) and i,j. Introducing the bounds 


A.. . = GLB., A (z ,<p) 

ij,mm 13 n 


A.. = LUB.. A (z , <p) 

13, max 13 n' 


4.. = 7-LUB,. A(z ,<p) , 4.. . = f GLB.. A(z,<p) (4.9) 

T3,max A 13 ' ^ 13, mm A 13 

G G 


p.. . = GLB.. p . (z ,03) , p.. = LUB.. p (z,<p) 

n.j,min 13 ^rnin ^ 13, max ij max ^ 
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taken over block i , (j — 1) and its neighbor i,j (as done for R^ in Section III), 
we can write 


A 3 £. n+2 . A.. . p.. . = R... . <R (z,<p)<R.. 

e i],mm ij,min 'ij.min ij,mm n' ij, max 


<a 3 e. n+2 

e ij , max 


A.. 

ij, max 


p.. 

'ij , max 


( 4 . 10 ) 


Estimates of . and can probably be made with reasonable ac- 

ij,min ij, max 

curacy. If it turns out that A is small, the estimation of n and 

ij , max "ij , max 

a. . can be based on knowledge of the density function near the surface of the 
ij , mm & 

Earth where considerable direct observational data is available. 


In any case, we see that A^(z , <p) rather than £^(z ,<p) itself is the sig- 
nificant parameter, and "block'' bounds, rather than point values, are what we 

need. As a first step, to get a feel for the order of magnitude of this parameter, 

7 

we might consider a recent radial model, constructed by Wang, for the density 

distribution of the Earth. In this model, the Earth is considered to be a sphere, 

the density depends only upon r , and a great deal of averaging is implicit in the 

construction. The basic idea would be to use Wang's model to approximate the 

two integrals of Eq. (4. 6): the first for various values of n, and the second for 

various values of the lower limit ( £ =1 for this model) , and then by direct 

max _ _ 

comparison, construct tables of values for £ n and A n =l-£ n . Following this, 
one should carry out similar calculations based on more realistic models and/or 
construct density variations along selected "rays" (z ,<p), for which a good supply 
of density data is available. Clearly, considerable analysis and numerical ex- 
perimentation will be required to determine realistic estimates of A.. 

ij, max 

and A.. 

ij, min 

It is possible that the theory of isostatic compensation can be used in the 
analysis. For if £ n (z,<p) turns out to be fairly close to, or below, the level of 
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compensation, it will follow that the second integral in Eq. (4. 6) will be quite in- 


sensitive to z and (o , so that the block bounds R.. . and R.. 

ij.min ij,max 

determined primarily by bounds for [ A(z , <p) ] . 


will be 
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V. "WORST" CASE ANALYSIS AND "BEST" RESULTS 


In the estimation of bounds by rigorous procedures, such as those imple- 
mented by Cholshevnikov, inequality relationships are used and the question arises 
as to whether the final result "really" contains a < or a ^ condition. In some 
cases, the method of analysis clearly indicates that equality is not a possibility, 
as in Cholshevnikov' s first and third theorems. If equality is not a possibility, 
then, clearly, the result is not "best": a sharper analysis should exist for which 
the bound can be tightened to the point where equality is a possibility. The realiza- 
tion of the condition of equality is a "worst case" because it, in effect, limits the 
further tightening of the bound. 

One way to test whether or not a result is "best" is to try to construct the 

"worst case" consistent with the hypotheses of the analysis. This is the method 

employed by Cholshevnikov in showing that his Theorem n is a "best" result for 

m = 0 . His "worst case" is a homogeneous hemisphere bounded by the z = 0 

3 

plane ("equatorial" plane) for which C decreases as l/n . This very extreme 
density distribution is consistent with the hypothesis of bounded density with bounded 
variation of his theorem. Actually, the numerical factors for C and the bounds 
of Theorem n do not quite match and neither does the exact dependence upon n ; 
what Cholshevnikov is really claiming is a "qualitatively best" result. 

Now, of course, we know that the Earth is not a homogeneous hemisphere; 

-3 

after all, the facts that J ~10 jj, and all other coefficients appear to be of order 
_6 “ 

±10 fj, indicate that the Earth is nearly a homogeneous sphere. If constraints are 
imposed on the density function so that a "worst case" construction is limited to 
more or less realistic density models, it is highly unlikely that any of the theorems 
quoted in Section II would qualify as best results. 
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One interesting conjecture is confirmed both by Cholshevnikov's example and 
by the less extreme construction outlined below. The conjecture is that a worst case 
constructed to match a bound, or more accurately to approximate a bound, for some 
subset of coefficients, will most likely imply that many other coefficients are far 
below their theoretical bounds. In Cholshevnikov's example, all the tesseral co- 
efficients C , S (m>0) vanish. A less extreme construction assumes realistic 
nm nm 

upper and lower bounds p . and p for the density. Then a particular C — 

^ 'min 'max nm 

(or S ) is maximized by assigning density p to all blocks for which the cor- 

nm max 

responding spherical harmonic is positive, and P min to the rest. It is easy to show 

that this distribution implies that S — (or C — ) vanishes, and that (at least) all 

nm ' nm' . 

C (or S ) with m - 2kin also vanish for k any integer > 0 . Similar results 
nm v nm 7 

can be obtained when an average value for p is imposed as an additional constraint. 

"Worst case" constructions are useful in a number of ways although, as in- 
dicated above, some care must be exercised in the interpretation of the implications. 
Certainly such constructions are useful to obtain a "feel" for the problem in the early 
stages of analysis, and it is in this application that caution must be observed since 
such constructions are usually based on a very few "global" parameters. In later 
stages of the analysis when, perhaps, some results have been obtained using a limited 
data base, construction of a worst case can be revealing: the extent to which the 
worst case satisfies constraints imposed by data, which are not incorporated in the 
analysis, could be a measure of the practical value of the results. Finally, a worst 
case construction utilizing an extensive data base is a possible means of obtaining 
realistic bounds, and, in fact, the methods outlined in Sections III and IV might be 
regarded as algorithms for such constructions. 



* 1 
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